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(57) Abstract: A method of iwninvasively 
determining die pulmonary capillary 
blood flow or cardiac output of a patient, 
including the use of rebreaming techniques 
to measure the respiratory flow and 
carbon dioxide pressure of the patient's 
breathing. These measurements are used 
to calculate carbon dioxide elimination 
and an indicator of the carbon dioxide 
content of the patient's blood, such as 
the end tidal pressure of carbon dioxide 
(PetCC>2), CaCOi, or pOOj. The location 
and orientation of a best-fit line through 
the carbon dioxide elimination data 
and the data of the indicator of carbon 
dioxide content is deternuned by linear 
regression or by plotting the carbon 
dioxide elimination data against the data 
of the indicator of carbon dioxide content 
At least one set of the data is modified 
and at least one other determination of 
the best-fit line made to find a data set that 
correlates most closely to the best fit line. 
The data may be modified by filtering or 
clustering. The slope of the best-fit line 
that correlates most closely to the data 
is then used to determine the pulmonary 
capillary blood flow or cardiac output of 
the patient When the indicator of carbon 
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isCaCOi. the negative of the slope of the best fit line is equal to the pulmonary capillary Wood flow of the patient. 
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NONINVASIVE DETERMINATION OP CARDIAC OUTPUT, PULMONARY BLOOD FLOW, AND BLOOD 
GAS CONTENT 

5 

TECHNICAL FIELD 
The present invention relates to methods and apparatus for accurately, 
noninvasively measuring the pulmonary capillary blood flow, ca r d i ac output, and mixed 
10 venous carbon dioxide content of the blood of a patient Particularly, the present 
invention relates to a method and apparatus for noninvasively measuring pulmonary 
capillary blood flow or cardiac output that employs an algorithm to increase the 
accuracy of data upon which the pulmonary capillary blood flow or cardiac output 
measurement is based. 

BACKGROUND 

Carbon dioxide elimination (VCQj) is the volume of carbon dioxide (COj) 
excreted from the body of a patient during respiration. Conventionally, carbon dioxide 
elimination has been employed as an indicator of metabolic activity. Carbon dioxide 
20 elimination has also been used in re-breathing methods of determining pulmonary 
capillary blood flow and cardiac output ' 
The carbon dioxide Fick equation: 

Q = Vcpj/CCvCOj - CaCOa), (1) 
where Q is cardiac output, CvC0 2 is carbon dioxide content of the venous blood of the 
25 patient, and CaC0 2 is the carbon dioxide content of the arterial blood of the patient, has 
been employed to non-invasively determine the pulmonary capillary blood flow or 
cardiac output of a patient The carbon dioxide elimination of the patient may be non- 
invasively measured as the difference per breath between the volume of carbon dioxide 
inhaled during inspiration and the volume of carbon dioxide exhaled during expiration, 
30 and is typically calculated as the integral of the carbon dioxide signal* or the fraction of - 
respiratory gases that comprises carbon dioxide, or "carbon dioxide fraction/* times the 
rate of flow over an entire breath. 

The partial pressure of end tidal carbon dioxide (PetC0 2 or etCOa) is also 
measured in rebreathing processes. The partial pressure of end-tidal carbon dioxide, 
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after correcting for any deadspace, is typically assumed to be approximately equal to the 
partial pressure of carbon dioxide in the alveoli (PaCOj) of the patient or, if mere is no 
intrapulmonary shunt, the partial pressure of carbon dioxide in the arterial blood of the 
patient (PaCQj). 

5 Rebreathing is typically employed either to non-invasively estimate the carbon 

dioxide content of mixed venous blood (as in total rebreathing) or to obviate the need to 
know the carbon dioxide content of the mixed venous blood (by partial rebreathing). 
Rebreathing processes typically include the inhalation of a gas mixture mat includes 
carbon dioxide. During rebreathing, the carbon dioxide elirnination of the patient 

10 decreases to a level less than during normal breathing. Rebreathing during which the 
carbon dioxide elirnination decreases to near zero is typically referred to as total 
rebreathing. Rebreathing that causes some decrease, but not a total cessation of carbon 
dioxide elimination, is typically referred to as partial rebreatbhig. 

Rebreathing is typically conducted with a rebreathmg circuit, wmch causes a 

15 patient to inhale a gas mixture that includes carbon dioxide. FIG. 1 schematically 
illustrates an exemplary rebreathing circuit 50 that includes a tubular airway 52 that 
communicates air flow to and from the lungs of a patient Tubular airway 52 may be 
placed in communication with the trachea of the patient by known intubation processes, 
or by connection to a breathing mask positioned over the nose and/or mourn of the 

20 patient A flow meter 72, which is typically referred to as a pneumotachometer, and a 
carbon dioxide sensor 74, which is typically referred to as a capnometer, are disposed 
between tubular airway 52 and a length of hose 60, and are exposed to any air that 
flows through rebreathing circuit 50. Both ends of another length of hose, which is 
referred to as deadspace 70, communicate with hose 60. The two ends of deadspace 70 

25 are separated from one another by a two-way valve 68, wttchrnayte 

direct the flow of air through deadspace 70. Deadspace 70 may also include an 
expandable section 62. A Y-piece 58, disposed on hose 60 opposite flow meter 72 and 
carbon dioxide sensor 74, facilitates the connection of an inspiratory hose 54 and an 
expiratory hose 56 to rebreathing circuit 50 and the flow communication of the 

30 inspiratory hose 54 and expiratory hose 56 withhose60. During inhalation, gas flows 
into inspiratory hose 54 from the atmosphere or a ventilator <not shown). During 
normal breathing, valve 68 is positioned to prevent inhaled and exhaled air from 
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flowing through deadspace 70. During rebieathing, valve 68 is positioned to direct the 
flow of exhaled and inhaled gases through deadspace 70. 

The re-breathed air, which is inhaled from deadspace 70 during rebreathing, 
includes air that has been exhaled by the patient (Le., carbon dioxide-rich air). 
5 Drofagtotdrebieathii^subs^ 

expired during the previous breath. Thus, during total rebreathing, the partial pressure 
of end-tidal carbon dioxide (PetC0 2 or etCOj) is typically assumed to be equal to or 
closely related to the partial pressure of carbon dioxide in the arterial (PaCOj), venous 
(PvCO,), or alveolar (PaCOj) blood of the patient. Total rebreathing processes are 

1 0 based on the assumption that neither pulmonary capillary blood flow nor the content of 
carbon dioxide in the venous blood of the patient (CvCOj), changes substantially cruring 
the rebreathing process. The partial pressure of carbon dioxide in blood may be 
converted to the content of carbon dioxide in blood by means of a carbon dioxide 
dissociation curve, where the change in the carbon dioxide content of the blood (CvC0 2 

15 - CaCOj) is equal to the slope (s) of the carbon dioxide dissociation curve multiplied by 
the measured change in end tidal carbon dioxide (PetCOj) as effected by a change in 
effective ventilation, such as rebreathing. 

In partial rebreathing, the patient inhales a rnixture of "fresh" gases and gases 
exhaled during the previous bream. Thus, the patient does not inhale a volume of 

20 carbon dioxide as large as the volume of carbon dioxide mat would be inhaled during a 
total rebreathing process. Conventional partial rebreathing processes typically employ a 
differential form of the carbon dioxide Fick equation to detennine the pulmonary 
capillary blood flow or cardiac output of the patient, which do not require knowledge of 
the carbon dioxide content of the mixed venous blood. This differential form of the 

25 carbon dioxide Fick equation considers measurements of carbon dioxide elimina ti o n, 
CvCOj, and the content of carbon dioxide in the alveolar blood of the patient (CaCOj) 
during both normal breathing and the rebreathing process as follows: 

Qpcwbd = VCQib-VCPtd <2> 

(CvCOaB - CvC0 2D ) - (CACO aB - CaC0 2D ), 

30 where VC0i B and Vco 2D are the carbon dioxide production of the patient during before 

rebreathing and during the rebreathing process, respectively, and CvC0 2B and CvC0 2D 

are the content of CO, of the venous blood of the patient before rebreathing and during 
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the rebreathing process, respectively. In using the differential Fick equation to calculate 
cardiac output, CaC0 2B and CaC0 2D , the contents of C0 2 in the arterial blood of Hie 
patient before rebreathing and during rebreathing, respectively, are substituted in 
equation (2) for the CaC0 2 measurements. 
5 Again, with a carbon dioxide dissociation curve, the meas^^ 

used to determine the change in content of carbon dioxide in the blood before and 
during the rebreathing process. Accordingly, the following equation can be used to 
determine pulmonary capillary blood flow or cardiac output when partial rebreathing is 
conducted: 

10 Q-AVccVsAPetCOj. O 

Alternative differential Fick methods of measuring pulmonary capillary blood 
flow or cardiac output have also been employed. Such differential Fick methods 
typically include a brief change of PetC0 2 and Vop, in response to a change in effective 
ventilation. This brief change can be accomplished by adjusting the respiratory rate, 

15 inspiratory and/or expiratory times, or tidal volume. A brief change in effective 

ventilation may also be effected by adding COj, either directly or by rebreathing. An 
exemplary differential Fick method that has been employed, which is disclosed in 
Gedeon, A. et al. in 18 Med & Biol. Eng. & Comput. 411-418 (1980), employs a period 
of increased ventilation followed immediately by a period of decreased ventilation. 

20 The carbon dioxide elirnination of a patient is typically measured over the 

course of a breath by the following, or an equivalent, equation: 

VC0j = /^Vx/cojtf, (4) 
where V is the measured respiratory flow and/c^ is the substantially simultaneously 
detected carbon dioxide signal, or fraction of the respiratory gases that comprises 

25 carbon dioxide or "carbon dioxide fraction." 

Due to the measured respiratory constituents upon which Vcp, and PetCOj 
calculations are made, Voo, typically responds to rebreathing about one bream before 
PetCQj for the same breath. Accordingly, a Vcpj signal may lead a PetC0 2 signal by 
about one breath. Thus, at a particular point in time, the VcOj and PetC0 2 signals do 

30 not correspond to one another. As these values are often used to noninvasively 

deterrnine pulmonary capillary blood flow or cardiac output, the lack of correspondence 



WO 01/62148 PCT/US00/24044 

-5- 

between these values may lead to inaccuracies in the pulmonary capillary blood flow or 
cardiac output determination. 

In addition, measurements that are taken during spurious breams, or breaths 
which do not provide information relevant to pulmonary capillary blood flow or cardiac 
5 output, may act as noise that introduces inaccuracy into the noninvasive pulmonary 
capillary blood flow or cardiac output determination. 

When equation (4) is employed to calculate the carbon dioxide elimination of 
the patient from the respiratory flow and carbon dioxide fraction measurements over an 
entire breath, such miscorrelation or noise-induced inaccuracies in either the expiratory 
10 flow, the inspiratory flow, or both may cause inaccuracies in the carbon dioxide 
elimination determination or inconsistencies between carbon dioxide elimination 
detenninations. 

Accordingly, there is a need for a method of accurately, noninvasivery 
calculating pulmonary capillary blood flow and cardiac output 

15 

DISCLOSURE OF INVENTION 

The present invention includes a method and apparatus for noninvasivery 

measuring pulmonary capillary blood flow and cardiac output The present invention 

includes the use of known rebreathing techniques to substantially noninvasivery obtain 

20 carbon dioxide eliiriination (VcOj) and partial pressure of end tidal carbon dioxide 

(PetCOj) measurements of a patient's breathing. These measurements may then be 

used to calculate pulmonary capillary blood flow or cardiac output of the patient by 

employing the following equation: 

Q -■ AVco, = AVCP? — (5) 
25 ACaC0 2 sAPetCO* 

where s is the slope of a standard carbon dioxide (COj) dissociation curve, AVcOj is the 

change in the carbon dioxide elimination of the patient due to a change in effective 

ventilation, such as that caused by rebreathing, and ACaCOj and APetC0 2 are the 

change in the content of carbon dioxide in the arterial blood of the patient and the 

30 change in the end tidal partial pressure of carbon dioxide of the patient, respectively, 

due to the same change in effective ventilation. Alternatively, a standard carbon 
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dioxide dissociation curve can be used to determine ACaCOj on the basis of the 
measured APetC0 2 . 

As an alternative to the use of the above equations to determine pulmonary 
capillary blood flow or cardiac output, the substantially noninvasive VcOj and CaCO a 
5 measurements can be related to each other in a linear fashion. This can be visually 
diagramed by plotting the Vcpj and CaC0 2 measurements against one another on a 
two-dimensional (X-Y) line graph. The negative slope (-1 * m) of the best-fit line 
through the data is approximately equal to the pulmonary capillary blood flow. The 
appropriate location and orientation of such a best-fit line maybe calcula^ by linear 

10 regression or least squares. Depending on the correlation between the calculated best- 
fit 1^ and the measured data, it may also be desirable to modify the data to provide a 
best-fit line that closely corresponds to Ihe data. 

In one embodiment of the method and apparatus of the present invention, the 
data can be modified by use of a known filter, such as a low-pass filter or a high-pass 

15 filter. Either digital or analog filters may be used. Eimer lhiear or nonlinear (e.g., 
median) filters maybe used. By way of example, and not to limit the scope of the 
present invention, a low-pass filter may be applied to the measured VcOj signal As 
another example, a high-pass filter may be applied to the measured CaC0 2 signal. 
Preferably, the filter and filter coefficient that are selected maximiz e die correlation 

20 between the measured VoOj and CaC0 2 signals. 

hi another embodiment of the method and apparatus of the present invention, 
the data points can modified by clustering. That is, the data points that are grouped 
closest to other data points are assumed to most accurately represent the true VcOj and 
CaC0 2 of the patient For example, the measured data with at least a predetermined 

25 number of close, or similar (e.g., within a specified threshold), data points is retained, 
while measured data with less man the predetermined number of close data points is 
discarded. The retained data points are assumed to be located on or near the best-fit 
line. In clustering, only these closely grouped sets of data points are considered in 
recalculating the best-fit line for the data and, thus, the negative slope (U^ -1 x m) of 

30 the best-fit line to determine the pulmonary capillary blood flow or cardiac output of 
the patient 
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Another embodiment of the method and apparatus of the present invention 
includes modifying the data points that are most likely to be closest to an accurately 
placed and oriented best-fit line. Each data point, which has a carbon dioxide 
elimination component (e.g., a y-ordinate component) and a component based on an 
5 indicator of carbon dioxide content (e.g., an x-ordinate component), is evaluated on the 
basis of a predetermined rninimum expected pulmonary capillary blood flow and a 
predetermined maximum pulmonary capillary blood flow. Lines, or the equations 
therefor, for both rninimum expected and maximum expected pulmonary capillary 
blood flows are located so as to intersect at each data point Then, the number of the 

10 other data points that are located between the two pulmonary capillary blood flow lines 
or equations is determined for each data point Only those data points with a threshold 
number of other data points between the two intersecting lines are used in the 
determination of the location and orientation of Ihe best-fit line through the data. 
Of course, any combination of methods of modifying data may be used to 

15 accurately determine the slope of the best-fit line through the measured VcOj and 
PetCOj data and, thus, to determine the pulmonary capillary blood flow or cardiac 
output of a patient 

The best-fit line through carbon dioxide elimination and carbon dioxide content 
data may also be used to determine the mixed venous carbon dioxide content of the 

20 patient when partial rebreaming techniques are employed to obtam flie data. As the 

mixed venous carbon dioxide content is assumed to equal the carbon dioxide content of 
the patient's blood when carbon dioxide elimination ceases (which does not occur 
during partial rebreathing), a best-fit line obtained by use of partial rebreatiring 
techniques can be used to noninvasively determine carbon dioxide content and, thus 

25 mixed venous carbon dioxide content, when carbon dioxide elimination is set at zero. 

Other features and advantages of the present invention will become apparent to 
those of ordinary skill in the art through a consideration of the ensuing description, the 
accompanying drawings, and the appended claims. 

30 BRIEF DESCRIPTION OF DRAWINGS 

FIG. 1 is a schematic representation of an exemplary rebreaming circuit mat 
may be employed with the methods of the present invention; 
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FIG. 2 is a schematic representation which illustrates the componentry that may 
be utilized to measure respiratory profile parameters that are employed mlhe methods 

of me present invention; 

FIG. 3 A illustrates an idealized, bi-directional rebreathing cycle with VoOj 
5 values for different breaths depicted as diamonds and PetCO a values for various breams 
shown as squares; 

FIG. 3B is a two-dimensional plot fflustrating the use of a known, H-<iirectional 
rebreathing process to obtain three Vcq, values and three values representative of the 
carbon dioxide content of the blood of a patient before, during, and after rebreathing; 
10 these three values have been used to substantially noninvasively determine the 
pulmonary capillary blood flow or cardiac output of the patient; 

FIG. 3C is a two-dimensional plot of a number of Vc^ values against the same 
number of carbon dioxide content values obtained over a single bi-directional 
rebreathing cycle; 

15 FIG. 3D is an exemplary two-dimensional plot depicting Vco 2 and carbon 

dioxide content values from the same re-breathing cycle as that shown in FIG. 3C and 
modified in accordance with the method of present invention; 

FIGs. 4A and 4B are two-dimensional plots mustrating an embodiment of a 
method for modifying data to obtain an accurate best-fit line therethrough in accordance 
20 with teachings of the present invention; 

FIG. 5 is a two-dimensional line graph mustrating a typical plot of VoOj an 
they-axis and CaC0 2 on the x-axis; and 

FIG. 6 is a two-dimensional line graph iUustrating a plot of Vop2 on the y-axis 
and CaC0 2 on the x-axis after the Vco, and CaC0 2 data have been modified in 
25 accordance with teachings of tiie jireserit invention. 

BEST MODES FOR CARRYING OUT THE INVENTION 
The present invention includes use of the Fick equation to calculate pulmonary 
capillary blood flow or cardiac output as the ratio of a change in carbon dioxide 
30 elimination, or VoOj, to a change in the content of carbon dioxide, or CaCO* in the 
arterial blood of a patient: 



WO 01/62148 



PCT/USOO/24044 



Q ° AVOO, (6) 
ACaCQj. 

CaCOj, or the contort of carbon dioxide in the arterial blood of a patient, can be 
noninvasively estimated by determining the PetCO* or partial pressure of carbon 
5 dioxide in the end tidal respiration of a patient and converting PetCOj to CaC0 2 by use 
of a standard carbon dioxide dissociation curve, as is known in the art, as follows: 

ACaC0 2 = sAPetCO„ CO 
where s is the slope of the carbon dioxide dissociation curve and APetC0 2 is a change 
in the end tidal partial pressure of carbon dioxide of a patient effected by a change in 
10 ventilation. Thus, pulmonary capillary blood flow or cardiac output can also be 
calculated as follows: 

Q = AVco 2 /sAPetCQ 2 . ( 8 ) 
Other indicators of the carbon dioxide content in the blood of a patient, such as 
pCOj, may be used in place of PetC0 2 or CaCOj to determine the pulmonary capillary 
15 blood flow or cardiac output of a patient 

VcOj and PetCOj, CaCO^, pCO^ or other indicators of the carbon dioxide 
content in the blood of a patient can be calculated or determined on the basis of 
substantially noninvasively obtained respiratory flow and respiratory carbon dioxide 
pressure data . 

20 HG. 2 schematically illustrates an exemplary method of substantially 

noninvasively monitoring the respiration of a patient and of measuring the flow rates 
and carbon dioxide concentration of gas mixtures that are inhaled and exhaled by a 
patient 10 over the course of the patient's breathing, such as during normal respiration 
or during known rebreathing techniques. A flow sensor 12 of a known type, such as the 

25 dirKawndal-pressme^peiespiratoiy 

Medical Systems Inc. ("Novametrix") of Wallingford, Connecticut (e.g., the 
Pediatric/Adult How Sensor (Catalog No. 6717) or the Neonatal Flow Sensor (Catalog 
No. 671 8)), which may be operatively attached to a ventilation apparatus (not shown), 
as well as respiratory flow sensors based on other operating principles and 

30 manufactured or marketed by others, may be employed to measure the flow rates of the 
breathing of patient 10. 
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A carbon dioxide sensor 14, such as the CAFNOSTAT® carbon dioxide sensor 
and a complementary airway adapter (e.g., the Pediatric/Adult Single Patient Use 
Airway Adapter (Catalog No. 6063), the Pediatric/Adult Reusable Airway Adapter 
(Catalog No. 7007), or the Neonatal/Pediatric Reusable Airway Adapter (Catalog No. 
5 7053)), which are manufactured by Novametrix, as well as main stream and side stream 
carbon dioxide sensors manufactured or marketed by others, may be employed to 
measure the carbon dioxide concentration of gas mixtures that are inhaled and exhaled 
by patient 10. 

Flow sensor 12 and carbon dioxide sensor 14 are connected to a flow monitor 
16 and a carbon dioxide monitor 18, respectively, each of which may be operatively 
associated with a computer 20 so that data from the flow and carbon dioxide monitors 
16 and 18 representative of the signals from each of flow sensor 12 and carbon dioxide 
sensor 14 may be detected by computer 20 and processed according to progr amming 
(e.g., by software) thereof. Preferably, raw flow and carbon dioxide signals from the 
flow monitor and carbon dioxide sensor are filtered to remove any significant artifacts. 
As respiratory flow and carbon dioxide pressure measurements are made, the 
respiratory flow and carbon dioxide pressure data may be stored by computer 20. 

Each breath, or breathing cycle, of patient 10 may be delineated as known in the 
art, such as by continuously monitoring the flow rate of the breathing of patient 10. 

As use of the Fick equation to calculate pulmonary capillary blood flow or 
cardiac output requires that a change in VcOj and CaCOj, PetCOj, pC0 2 or another 
indicator of the carbon dioxide content in the blood of a patient be known, a change in 
effective ventilation is required. By way of example, and not to limit the scope of the 
present invention, rebreathing techniques, such as by use of a dead space 70 such as that 
provided by the rebreathing circuit illustrated in FIG. 1, may be employed to cause a 
change in effective ventilation. FIG. 3 A illustrates the changes that may occur when a 
bi-directional rebreathing process is used to effect a change in effective ventilation. 
The graph of FIG. 3A illustrates the typical changes in the VcOj (shown as diamonds) 
and carbon dioxide content measurements (e.g., PetCO* shown as squares) that may 
occur between the baseline breathing (Le,, before rebreathing), during rebreathing, and 
the stabilization (le., after rebieathing) periods of an idealized {Le., without noise) 
bi-directional rebreathing cycle. During rebreathing, Vopj changes from a baseline 
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value (e.g., about 200 ml/min) to a during rebreathing plateau (e.g., of about 100 
ml/min.) within about 3 or 4 breaths, whereas carbon dioxide content may take longer 
to change from a baseline value (e.g., 38 mmHg) to a plateau (e.g., about 35 mmHg). 
FIG. 3B is a two-dimensional plot illustrates that one value, the plateau value, 
5 from each of the before, during, and after rebreathing phases of a bi-directional 
rebreathing process, such as that illustrated in FIG. 3 A, were used to estimate 
pulmonary capillary blood flow or cardiac output By way of contrast, in a method of 
determining pulmonary capillary blood flow or cardiac output incorporating teachings 
of the present invention, VcOj and carbon dioxide content data are continually 
1 0 measured, providing a plot such as that shown in FIG. 3C, with data at 1 00 being based 
on before rebreathing measurements, data along arrow 102 being based on during 
rebreathing measurements, and data along arrow 104 being based on after rebr eathing 
measurements. These data may be obtained by use of a single rebreathing cycle, over 
the course of a number of rebreathing cycles, at one or more discrete time intervals, or 
15 on a breath by breath basis, where data is continually measured, calculated, and 

analyzed in accordance with the method of the invention so as to continually update or 
monitor the pulmonary capillary blood flow or cardiac output of a patient 

When rebreathing or other known techniques are used to cause a change in 
effective ventilation so as to facilitate the substantially noninvasive detennination of 
20 pulmonary capillary blood flow or cardiac output, respiratory flow and carbon dioxide 
pressure data are obtained during at least the before, during, and after stages of 
rebreathing. Total or partial rebreathing processes may be used in the method of the 
present invention. These respiratory flow and carbon dioxide pressure data are then 
used, as known in the art, to calculate Vco 2 and PetCOj, as well as the changes in VCQj 
25 and PetC0 2 that occur with the change in effective ventilation. 

The calculated VcOj and PetC0 2 data are then used to determine the pulmonary 
capillary blood flow or cardiac output of the patient, such as by use of the Fick 
equations presented above. 

As an alternative, the pulmonary capillary blood flow or cardiac output of a 
30 patient can be determined over the course of a plurality of breaths by expressing the 
calculated Vcd^ data and CaC0 2 data or data of another indicator of the content of 
carbon dioxide in the blood of a patient, such as PetC0 2 or pCOj, in two-dimensions, 
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such as on a two-dimensional (X-Y) line graph, with Vcpj data points being measured 
on the y-axis and PetCO a data points being measured on the x-axis, then identifying a 
line that best fits the data, which is also referred to herein as a best-fit line. - 
For example, me equation for the best fit line is: 
5 y = mx + b (9) 

or 

m-v-b 00) 
x , 

where y is the y-axis coordinate of a data point, x is the x-axis coordinate of the same 

10 data point, m is the slope of the line, and bis the ofiset value for the line. If VoOj is 

measured on the y-axis and CaC0 2 is measured on the x-axis, then 

m^ VCQ 7 -b (ID 
CaCOj . 

The negative slope (i.e., -1 * m) of the best-fit line through the Vcc^-CaCC^ data would 
15 be equal to the pulmonary capillary blood flow or cardiac output of the patient: 

-m = Q. d2) 
The best-fit line for the VcOj and CaCO z data is preferably deteimined by use of 
known linear regression techniques or any other known methodology for determining 
the relationship between two variables. The method of linear regression provides an 
20 accurate pulmonary capillary blood flow or cardiac output value based on a large 

number of VcOj and CaC0 2 data obtained over the course of one or more changes in 
effective ventilation. When linear regression is used, the slope (m) of the best-fit line 
for the data is calculated as follows: 

m = Lxy/Lxx (I 3 ) 
25 and the offset (b) of the line is calculated by the following equation: 

b = Ey/n — m * Exm, (14) 

where 

Lxx = 2x a -(Exx2x)/rv (I 5 ) 
Lyy=2Jy 2 _(^y xSyyn,and (1$) 
30 Lxy=Exy -<Ex xEyVn, (17) 

and where n is the number of data points in the plot, Ex is the sum of all x-<xx)rdinate 
(LC CaCO a content) values, Ey is the sum of all y-coordinate Ci.e., VcOj) values, Ex 2 is 
the sum of the square of all xK»ordinate values, Ey 2 is the sum of me square of all y- 
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coordinate values, and Sxy is the sum all paired x- and y-coordinate values multiplied 
by each other. 

When linear regression is used to determine the location and orientation of a 
best-fit line, a correlation coefficient (r) that quantifies the accuracy -with which the 
5 best-fit line correlates to the VcOj and CaC0 2 data can also be calculated as follows: 

r = (Lxy x Lxy)/(Lyy * Lxx). (18) 
Alternatively, any other measure of the quality of fit that quantifies the accuracy with 
which the best-fit line correlates to the VCOj and CaC0 2 data may be used- 
Correlation coefficients range from 0 to 1 .0, where a correlation coefficient of 0 
10 indicates that no linear correlation exists between the x-ordinate and the y-ordinate data 
and a correlation coefficient of 1 .0 indicates that the x-ordinate and y-ordinate data are 
perfectly linearly correlated (i.e., all of the Vcx^-CaCOj data points are located on the 
same straight line). 

The Vcc^-CaCOj data points measured before and during rebreathing, however, 
15 are rarely located on the same straight line. One reason forinis is that, during 

rebreathing maneuvers, the Vopj signal typically leads the PetC0 2 signal and, thus, the 
CaC0 2 by about one breath. In addition, VcOj is calculated on the basis of signal 
components that have higher frequencies than do the PetC0 2 signal. As a result, when 
the VcOj and CaCO a measurements calculated over a period of time are plotted against 
20 one another on a two-dimensional (X-Y) line graph, the result typically appears as an 
arc or a loop, as shown in FIGs. 3C and 5, rather than as a straight line, depending on 
the amount of data calculated and the duration of rebr eathing . Moreover, VcOj and 
CaC0 2 measurements may be calculated on the basis of respiratory flow and carbon 
dioxide pressure data obtained during spurious breaths. Such data do not relate to the 
25 pulmonary capillary blood flow or cardiac output measurement. Vc^andCaCOj 
calculations that are based upon such spurious date 

miscalculation of a best-fit line through the calculated VcOj and CaC0 2 data. Asa 
result, the correlation coefficient of a best-fit line to the data is typically much less than 
1.0. 

30 The measured respiratory flow and carbon dioxide pressure data or die 

calculated Vcpj and CaC0 2 data can be modified to increase the correlation coefficient 
between the Vco 2 and CaC0 2 data and the best-fit line therefor. Preferably, a linear 
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transform is used to increase the correlation coefficient A linear transform may be 
used to delay the calculation of a Vco 2 data point to accurately coincide therewith a 
CaC0 2 data point based on measurements taken during the same bream. The measured 
or calculated data may also be filtered by use of a linear transform. 
5 In one embodiment of a method for increasing the correlation coefficient 

between the VcOj and CaC0 2 data and the best-fit line therefor, a filter is applied to the 
calculated VoOj or CaCO a data. Known analog or digital low-pass, high-pass, or band 
pass filters, including adaptive filters, may be employed. Linear or nonlinear may be 
employed. Preferably, a first order (single pole) infinite impulse response (HR) digital 
10 filter is employed to filter the VCQj calculations in a manner that improves the 

correlation between the Vco 2 calculation and the lagging PetC0yCaCO 2 calculation. 

The equation for such a filter is: 

Vco 2 , [n] = oxVcp 2 < [n-l] + (l-«)'<Vco 2 [n], (19) 
where Vcqjn] is the most recently calculated, unfihered Vcdj data point, Vc<V£n-l] is 

15 the previous, filtered VcCfe data point, VcOj' [n] is the new "filtered" value based on 
VCOjfn] and obtained by use of the filter, and a is the filter coefficient The filter 
coefficient, a, has a range of 0 to 1.0. The greater the value of o, the more profoundly 
the most recently calculated data point is filtered and, conversely, the lower a values 
cause the most recently calculated data points to be filtered to a lesser degree. When a 

20 is equal to zero, the most recently calculated data point is not filtered. 

Due to anatomical and physiological differences between different patients, 
different patients have differing optirnal filter coefficients, a. In addition, as anatomical 
and physiological changes may occur in a patient over time, the optimum filter 
coefficients, o, to be used in filtering the VcOj or CaC0 2 values calculated from the 

25 patient's breathing may also vary over time. Accordingly, the selection of an optimal 
filter coefficient, a, is also within the scope of the present invention. Any known 
optimization method or search algorithm may be employed to select optirnal filter 
coefficient, a. 

As an example of the way in which an optimal filter coefficient may be selected, 
30 a is first set to a default value (e.g., 0.85) and the calculated VcOj or,CaC0 2 values are 
filtered on the basis of the default filter coefficient, a. The linear regression is then 
performed to obtain a best-fit line. If the correlation coefficient of best-fit line 
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calculated with the just-filtered data is less than the correlation coefficient of the 
immediately preceding best-fit line, which was calculated with unfiltered data or with a 
prior filter coefficient, then a predetermined a adjustment value (e.g., 0.01) is changed 
by multiplying the a adjustment value by -1 and by modifying the filter coefficient by 
5 adding the modified a adjustment value thereto. Otherwise, the filter coefficient, a, is 
modified by adding the unmodified a adjustment value thereto. The process of filtering 
the data based on a modified filter coefficient, obtaining a best-fit line for the data, 
comparing the correlation coefficient of the best-fit line to the correlation coefficient of 
the previous best-fit line, and adjusting the filter coefficient accordingly is then repeated 
1 0 a predetennined number of times (e.g., 50 times). The best-fit line with the greatest 
correlation coefficient, based on the unfiltered data and each set of filtered data, is 
selected to calculate the pulmonary capillary blood flow or cardiac output of the patient 
When filtering is used, the Vc0i-CaCO 2 plot preferably narrows, as depicted in 
FIGs. 3D and 6, to thereby increase the accuracy with which fee location and 

15 orientation of a best-fit line can be established and, thus, to increase the accuracy of a 
pulmonary capillary blood flow or cardiac output determination based on the data. 

Another embodiment of a method for increasing the correlation coefficient 
between the Vcpj and CaC0 2 data and the best-fit line therefor, which is referred to 
herein as "clustering," includes the selection of data points that are grouped closely 

20 together. That is, the data points that are selected include those data points having a 
number of other data points within a predetennined range thereof Data points that are 
not clustered are probably inaccurate or based on measurements taken during spurious 
breaths. As an accurate best-fit line through the data would likely be based on the 
clustered data, the data points that are not located in a cluster are not used in calculating 

25 the location and orientation of a best fit line for the data. 

Clustering of the data points may include normalization or transformation of the 
data such that ranges of the x-coordinate data (e.g., the CaC0 2 data) and the 
y-coordinate data (e.g., the VcOj data) are substantially the same. Without such 
normalization, the data group (e.g., the VcOj data or the CaC0 2 data) with the highest 

30 range would dominate; the other data group would be less significant. 

An exemplary manner in which the data may be normalized includes use of the 
following normalization: 
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x = (x-5?yo x , (20) 

where: 

x is the raw value, 5? is the mean value of all x-axis (e.g., CaCOj) data in the plot, and 
a x is the standard deviation of all x-axis data in the plot This normalization is applied 
5 to all x-axis values. A similar normalization scheme is applied to all of the y-axis 
values. 

The normalized data may then be clustered by searching for a redetermined 
number (e.g., 5) of the closest data points (e.g., VcOi or CaC0 2 data points) to each of 
the data points in a group. The differences between the analyzed data point and each of 

10 the predetermined number of closest data points are men added together and compared 
to a predetermined threshold. If the sum of the differences exceeds the predetermined 
threshold, the analyzed data point is discarded. Of course, the use of other clustering 
techniques to identify the most accurate data and to disregard probable inaccurate data 
are also within the scope of the present invention. 

15 Once clustering has been performed, the inverse of the no rmaliza ti o n is 

calculated, or the normalization is undone, to provide an accurate determination of 
pulmonary capillary blood flow or cardiac output An example of the manner in which 
the inverse of the normalization may be calculated includes use of the following 
equation: 

20 x = xa»+5t (21) 

This inverse of the normalization is applied to all of the clustered x-axis (e.g., CaCOj) 
values. A similar inverse normalization scheme is applied to all of the clustered y-axis 
date. 

Clustering is one of many known techniques for detennini^ Other 
25 known techniques for determining outliers may also be used in the method of the 
present invention. 

Alternatively, or in addition to disregarding probable inaccurate data points, in 
order to enhance the accuracy of the data, clustering can be used add synthetic data 
points. Synthetic data points may be added to increase the correlation coefficient of Ac 
30 best-fit line to the data points on which the best-fit line is based. 

Another embodiment of the method for modifying data that incorporates 
teachings of the present invention is depicted in FIGs. 4A and 4B. As with the filtering 
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and clustering embodiments described previously herein, the present embodiment 
includes selection of data points that are most likely to facilitate an accurate 
determination of the location and orientation of a best-fit line and, thus, of the 
pulmonary capillary blood flow or cardiac output of a patient. This embodiment of the 
5 method for modifying data includes iteratively examining data points and the 

distribution of the remaining data points relative to the two lines representing the range 
of possible PCBF measurements. 

As shown in FIGs. 4A and 4B, a line or the equation for a line 1 10 representing 
a mi™'™"™ expected pulmonary capillary blood flow (i.e., -10^= PCBF^ and a line 

10 or the equation for a line 120 representing a m a x i mum expected pulmonary capillary 
blood flow (i.e., -m llne m - PCBF^J are positioned to intersect at a data point 130. For 
example, when the x-ordinate is based on CaCO* line 1 10 may have a slope of -03, 
which represents a minimum expected pulmonary capillary blood flow of 0.5 LAnin, 
and line 120 may have a slope of -20, which represents a maximum pulmonary 

15 capillary blood flow of 20 L/min. Of course, the use of other pulmonary capillary blood 
flow values for lines 1 10 and 120 is also within the scope of the present invention. 

Next, the number of other data points 130 located between lines 1 10 and 120 is 
determined. If the number of data points 130 between lines 1 10 and 120 is equal to or 
exceeds a threshold number, the analyzed data point 130 is retained for a subsequent 

20 determination of the location and orientation of a best-fit line through the data. 

Otherwise, the analyzed data point 130 is discarded. The threshold number of data 
points that must be located between line 1 10 and line 120 for an analyzed data point to 
be retained may be a predetermined value or determined by other means. As example, 
the threshold number may be set to the median number of data points mat are located 

25 between line 1 10 and line 120 when each data point 130 of a set of data points 130 has 
been evaluated in accordance with the present embodiment of the method for modifying 
data. This process is repeated until each data point 130 in a set of data points 130 has 
been so evaluated. FIG. 4A depicts use of the present embodiment of the data 
modification method on a data point 130 that will be retained, while FIG. 4B illustrates 

30 use of the present embodiment of the data modification method on another data point 
130' that will not be retained. 
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FIGs. 3C and 3D and FIGs. 5 and 6 illustrate the effect of modifying data in 
accordance with teachings of the present invention to increase the accuracy with which 
the location and orientation of a best-fit line through the data may be determined. 
FIG. 5 illustrates a typical Vcpj vs. CaC0 2 plot without such modification, where the 
5 plot appears as a loop. By way of contrast, FIG. 6 illustrates the closeness of the data 
when one or more of the embodiments of the method of the present invention are used 
to modify the data. FIGs. 3C and 3D illustrate plots of VcOj and PetC0 2 data before 
and after modification in accordance with the present invention, respectively. The 
increased closeness of the data points makes it possible to determine the orientation and 
location of a best-fit line therethrough with increased accuracy. 

Once all of the data points have been examined, the location and orientation for 
the best-fit line through the remaining, clustered data are determined. Again, linear 
regression is preferably used to determine the location and orientation of the best-fit 
line. The negative slope (Le., -1 * m) of the best-fit line provides a pulmonary capillary 
blood flow measurement, which may then be used to determine cardiac output A 
correlation coefficient can then be calculated, as previously disclosed herein, to indicate 
the quality of the data used to determine pulmonary capillary blood flow or cardiac 
output The correlation coefficient or a quality measure based thereon may then be 
communicated to the user (e.g., a doctor, nurse, or respiratory technician) or used to 
weight the resulting pulmonary capillary blood flow or cardiac output value in an 
output weighted average value. 

One or a combination of the embodiments of the method for modifying data in 
accordance with the present invention may be performed on the measured or calculated 
data to increase the accuracy with which a best-fit line through the data or the 
pulmonary capillary blood flow or cardiac output of a patient can be det ermin ed. 

As an example of the use of filtering and clustering together, the calculated 
Vco 2 data are grouped together as the y-axis data of a two-dimensional line graph and 
the calculated CaC0 2 data points are grouped together as x-axis data points. The data 
points in at least one of the groups are filtered to determine a best-fit line for the data 
having an optimum correlation coefficient The data are also clustered, either before or 
after filtering, to further improve the correlation coefficient of the best-fit line to the 
calculated VcOj and CaC0 2 data. The remaining data is then used to determine (e.g., 
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by linear regression) a best-fit line therefor, as well as a correlation coefficient for the 
best fit line. The slope of the best fit line is then calculated and used to determine 
pulmonary capillary blood flow or cardiac output The correlation coefficient may also 
be used to indicate flie reliability of the pulmonary capillary blood flow or cardiac 
5 output detennination or to impart a specific weight to the pulmonary capillary blood 
flow or cardiac output determination in a weighted average thereof 

Once the location and orientation of an accurate best-fit line for the data has 
been determined, as disclosed previously herein, the pulmonary capillary blood flow of 
the patient can be calculated as the negative of the slope of the best fit line. In addition, 

10 cardiac output can then also be determined by adding the pulmonary capillary blood 
flow of the patient to the mtrapulmonary shunt flow of the patient, which can be 
determined by known processes. 

In addition, the best-fit line can be used to estimate mixed venous carbon 
dioxide content of the patient Conventionally, total rebreathing techniques have been 

1 5 required to substantially noninvasively measure mixed venous carbon dioxide content 
When carbon dioxide elirnination eventually ceases during total rebreathing, the partial 
pressure of carbon dioxide measured at the mouth of a patient may represent the mixed 
venous carbon dioxide content of the patient When partial rebreathing techniques are 
used, the carbon dioxide elimination of the patient is reduced to levels lower man 

20 baseline, but is not reduced to zero. By employing teachings of the present invention to 
determine the best-fit line through data obtained by use of partial rebreathing 
techniques, the best-fit line can be extended to a point where carbon dioxide 
elimination would be equal to zero or effectively zero andtoerebytodeterniinethe 
carbon dioxide content, or mixed venous carbon dioxide content, of the patient's blood 

25 at that point Equation (11), which is me equation for the best-fit line, can be 
rearranged in terms of carbon dioxide elimination as follows: 

VcOj = m x CaC0 2 + b. (22) 
When carbon dioxide elirnination ceases, VcOj is equal to zero and equation (22) 
becomes: 

30 0 = mx CvC0 2 +b, (23) 

where CvC0 2 is the mixed venous carbon dioxide content, which can be rearranged as 
follows: 
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CvC0 2 = -b/m. . (24) 

Accordingly, the present invention also includes a method for substantially 
noninvasively deteainining mixed venous carbon dioxide content when partial 
rebreathing techniques are employed 

5 

Although the foregoing description contains many specifics, these should not be 
construed as limiting the scope of the present invention, but merely as providing 
illustrations of some of the presently preferred embodiments. Similarly, other 
embodiments of the invention may be devised which do not depart from the spirit or 

10 scope of the present invention. Features from different embodiments may be employed 
in combination. The scope of the invention is, therefore, indicated and limited only by 
the appended claims and their legal equivalents, rather than by the foregoing 
description. All additions, deletions and modifications to the invention as disclosed 
herein which fall within the meaning and scope of the claims are to be embraced 

15 thereby. 
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CLAIMS 

What is claimed is : 

1. A method for noninvasively determining at least one of a mixed venous 
carbon dioxide content, a pulmonary capillary blood flow, and a cardiac output of a 

5 patient, comprising: 

determining a plurality of data comprising carbon dioxide e liminati on data and data of 

an indicator of carbon dioxide content in blood of the patient; 
determining a correlation coefficient between said carbon dioxide elimination data and 
said data of said indicator of carbon dioxide content; and 
10 calculating at least one of the mixed venous carbon dioxide content, the pulmonary 
capillary blood flow, and the cardiac output based on said correlation 
co efficient 

2. The method of claim 1, wherein said deter mining said correlation 

1 5 coefficient comprises plotting said carbon dioxide elimination data against said data of 
said indicator of carbon dioxide content 

3. The method of claim 2, wherein said plotting comprises plotting said 
carbon dioxide elimination data along a y-axis of a two-dimensional line graph against 

20 corresponding data of said indicator of carbon dioxide content plotted along an x-axis 
of said two-dimensional line graph. 

4. The method of claim 3, wherein said determining said correlation 
coefficient further comprises determining a best-fit line through the plotted plurality of 

25 data 

5. The method of claim 4, wherein said calculating the pulmonary capillary 
blood flow includes multiplying a slope of said best-fit line by -1 . 

30 6. The method of claim 1, further comprising modifying at least one of said 

carbon dioxide elimination data and said data of said indicator of carbon dioxide 
content to increase said correlation coefficient. 
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7. The method of claim 6, wherein said modifying comprises selecting data 
to use in said determining said correlation coefficient, said selecting comprising at least 
one of: 

filtering at least one of said carbon dioxide elimination data and said data of said 
5 indicator of carbon dioxide content; 

clustering said carbon dioxide elimination data and said data of said indicator of carbon 
dioxide content; and 

determining a number of data points between a first line representing a minimum 

expected mixed venous carbon dioxide content, pulmonary capillary blood flow, 
10 or cardiac output and a second line intersecting said first line at an analyzed data 

point of said plurality of data, said second line representing a maximum 
expected mixed venous carbon dioxide content, pulmonary capillary blood flow, 
or cardiac output and comparing said number to a threshold number in 
determining whether to retain said analyzed data point 

15 

8. The method of claim 7, wherein said filtering is effected on the basis of 
an infinite impulse response form employing the foUowing equation: 

Vcpj'[n] = o x VaV[n-l] + (l-<x) * VeOjfcn.], 
w here Vcojn] is the most recently obtained, unfiltered carbon dioxide elimination data 
20 point, VaV[n-l] is the previous, filtered carbon dioxide elimination data point, 

' Vco/M is the new "filtered" value based on Vco 2 [n] and obtained by use of the filter, 
and a is a filter coefficient. 

9. The method of claim 7, further comprising normalizing at least a portion 
25 of said plurality of data prior to said clustering. 

10. The method of claim 9, further comprising reversing said normalizing 
following said clustering and before said detemuning said correlation coefficient 



30 



11. The method of claim 1, wherein said deterniining said plurality of data 
comprising data of an indicator of carbon dioxide content in blood of the patient 
comprises determining CaC0 2 data or end tidal partial pressure data. 



wo 
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12. A respiratory evaluation apparatus for noninvasively determining at least 
one of a mixed venous carbon dioxide content, a pulmonaiy capillary blood flow, and a 
cardiac output of a patient, comprising: 

sensing means for obtaining a plurality of data comprising carbon dioxide elimination 
5 data and data of an indicator of carbon dioxide content in blood of the patient; 

first processing means for determining a correlation coefficient between said carbon 
dioxide elimination data and said data of said indicator of carbon dioxide 
content; and 

second processing means for calculating at least one of the mixed venous carbon 
10 dioxide content, the pulmonary capillary blood flow, and the cardiac output 

based on said correlation coefficient 

13. The respiratory evaluation apparatus of claim 12, wherein said first 
processing means is configured to determine said correlation coefficient by plotting said 

1 5 carbon dioxide elimination data against said data of said indicator of carbon dioxide 
content 

14. The respiratory evaluation apparatus of claim 13, wherein said first 
processing means is configured to plot said carbon dioxide elimination data along a y- 

20 axis of a two-dimensional line graph against corresponding CaC0 2 or end tidal partial 
pressure data plotted along an x-axis of said two-dimensional line graph. 

15. The respiratory evaluation apparatus of claim 14, wherein said first 
processing means is configured to determine said correlation coefficient by dete rminin g 

25 a best-fit line throiigh the plotted plurality of data. 

16. The respiratory evaluation apparatus of claim 15, wherein said second 
processing means is configured to effect said calculating the pulmonary capillary blood 
flow by multiplying a slope of said best-fit line by -L 



30 



17. The respiratory evaluation apparatus of claim 12, further comprising 
improvement means for modifying at least one of said carbon dioxide e liminatio n data 
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and said data of said indicator of carbon dioxide content to increase said correlation 
coefficient 

1 8. The respiratory evaluation apparatus of claim 17, wherein said 

5 improvement means is configured to select data to use in dete rminin g said correlation 
coefficient and includes at least one of: 

refinement means for filtering at least one of said carbon dioxide elimination data and 
said data of said indicator of carbon dioxide content; 

grouping means for clustering said carbon dioxide elimination data and said data of said 
10 indicator of carbon dioxide content; and 

evaluation means for determining a number of data points between a first line 

representing a minimum expected mixed venous carbon dioxide content, 
pulmonary capillary blood flow, or cardiac output and a second line intersecting 
said first line at an analyzed data point of said plurality of data, said second line 
j 5 representing a maximum expected mixed venous carbon dioxide content, 

pulmonary capillary blood flow, or cardiac output and comparing said number 
to a threshold number in detennining whether to retain said analyzed data point 

19. The respiratory evaluation apparatus of claim 18, wherein said 
20 refinement means is configured to effect said filtering on the basis of an infinite 

impulse response form by employing the following equation: 
VcOj^n] = a x VcOj^n-l] + (1-a) * Vcojn], 
where Vco 2 [n] is the most recently obtained, unfiltered carbon dioxide elimination data 
point, VaV[n-l] is the previous, filtered carbon dioxide elimination data point, 
25 VC(V[n] is new "filtered" value based on Vco 2 [n] and obtainedbyuseof the filter, 
and a is a filter coefficient 

20. The respiratory evaluation apparatus of claim 18, wherein said grouping 
means is further configured to normalize at least a portion of said plurality of data prior 

30 to said clustering. 
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21 . The respiratory evaluation apparatus of claim 20, wherein said grouping 
means is further configured to reverse said normalizing following said clustering and 
before said determining said correlation coefficient 

5 22. The respiratory evaluation apparatus of claim 12, wherein said data of an 

indicator of carbon dioxide content in blood of the patient comprises CaCO a data or end 
tidal partial pressure data. 



10 



23. 

processor. 



The respiratory evaluation apparatus of claim 12, further including a 
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